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ABSTRACT 

We show that evolution of the luminosity density of galaxies in the universe 
provides a powerful test for the geometry of the universe. Using reasonable 
galaxy evolution models of population synthesis which reproduce the colors of 
local galaxies of various morphological types, we have calculated the luminosity 
density of galaxies as a function of redshift z. Comparison of the result with 
recent measurements by the Canada-France Redshift Survey in three wavebands 
of 2800A, 4400A, and 1/im at z < 1 indicates that the A-dominated flat universe 
with A ~ 0.8 is favored, and the lower limit on A yields 0.37 (99% C.L.) or 0.53 
(95% C.L.) if Q + Ao = 1. The Einstein-de Sitter universe with (Oq, Ao) = (1, 0) 
and the low-density open universe with (0.2, 0) are however ruled out with 
99.86% C.L. and 98.6% C.L., respectively. The confidence levels quoted apply 
unless the standard assumptions on galaxy evolution are drastically violated. 

We have also calculated a global star formation rate in the universe to be 
compared with the observed rate beyond z ~ 2. We find from this comparison 
that spiral galaxies are formed from material accretion over an extended period 
of a few Gyrs, while elliptical galaxies are formed from initial star burst at z > 5 
supplying enough amount of metals and ionizing photons in the intergalactic 
medium. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: 
intergalactic medium — cosmology: theory 
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1. Introduction 

Recent progress in deep redshift surveys of galaxies reaching z ~ 1 gives valuable 
information on the history of galaxy evolution (Lilly et al. 1995; Ellis et al. 1996). 
The Canada-France Redshift Survey (CFRS) consisting of 730 /As-selected galaxies with 
17.5 < Iab < 22.5 revealed a marked evolution in their comoving luminosity density in 
three wavebands of 2800A, 4400A, and 1/iin in the range of z = 0-1 (Lilly et al. 1996). 
They argued that their data can be explained by a galaxy evolution model of population 
synthesis with star formation rate proportional to t~ 2 ' 5 where r is the elapsed time from the 
initial turn-on of star formation. In practice, however, their use of a single index —2.5 for 
all galaxies is not justified because the present-day colors of galaxies are strongly correlated 
with their morphological type through type- dependent history of star formation (Tinsley 
1980), and a single-type analysis is known to produce a bias in interpreting the data when 
applied to a sample with a realistic mixture of morphological types (Yoshii & Takahara 
1989). One therefore has to construct type- dependent evolution models of galaxies and sum 
up with relative proportion of their types. If the global history of star formation for the 
composite of all types of galaxies is determined in this way, cosmological models can then 
be tested against the observations of the luminosity density evolution in the universe. 

In this Letter, we calculate an evolution of the luminosity density using the galaxy 
evolution models of population synthesis developed by Arimoto & Yoshii (1987, hereafter 
AY) and Arimoto, Yoshii, & Takahara (1992, hereafter AYT) where the effect of chemical 
evolution is appropriately taken into account. Since the integrated color of stars in a galaxy 
is very sensitive to an assumed history of star formation, the present-day spectral energy 
distribution (SED) in the ultraviolet through near-infrared region for each type of galaxies 
can only be reproduced by almost identical luminosity evolution as a function of look-back 
time until the epoch corresponding to an average age of stars in a galaxy (AY, AYT). These 
models are therefore reliable especially at low redshifts of z < 1 (Yoshii & Takahara 1988) 
where we constrain the cosmological parameters from the CFRS data. A global rate of star 
formation is also calculated and is used to discuss an early evolution of galaxies from some 
recent observations of star formation rate at z > 1. 

2. Theoretical Model of the Luminosity Density Evolution 

Let £(t, A) be the comoving luminosity density of galaxies measured at time t from 
the epoch of galaxy formation and at rest-frame wavelength A per unit wavelength. Then 
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£(t, A) is expressed as, after summing over galaxy types, 

C(t, A) = E / dL B ML B )L k (t, A; L B ) , (1) 

k 

where ip k (LB) is the present-day luminosity function for blue-selected galaxies of /c-th type 
and L k (t, A; Lb) is the luminosity at t and A per unit wavelength for a galaxy of k-th type 
whose present-day 5-band luminosity is Lb- Assuming that the SED evolution f k (t,\) 
depends only on the galaxy type and not on the luminosity of each galaxy, we write 



L k (t, A; L B ) — 



fk(t,\) 
fk(T G , \ B ) 



L k (T G ,X B ;L B ) , (2) 



where T G is the present age of galaxies and A B =4400A. We note that the assumption of 
luminosity-independent evolution is consistent with the method of evaluating the luminosity 
density from the CFRS data where only galaxies brighter than ~ L* are sampled at z > 0.5 
and the luminosity function is extrapolated to fainter magnitudes with its shape kept 
constant (Lilly et al. 1996). Thereby, the luminosity density at t and A can be related to 
the present-day £>-band luminosity density of each type: 



£(i,A) = E 



k 



AM 

fk(T G , A B ) 



C k (T G , A B ) • (3) 



We classify galaxies into E/SO, Sab, Sbc, Scd, and Sdm which evolve with different star 
formation histories. The evolving SEDs of spiral galaxies are given in AYT over a range of 
wavelength spanning from 1550 to 34000A, and those for elliptical galaxies are taken from 
the 1O 12 M (baryon mass) model in AY for 3600-34000A. The ultraviolet SED of elliptical 
galaxies, the origin of which is still a matter of debate, is prescribed in the same way as in 
Yoshii & Peterson (1991), but the following results are hardly changed by this uncertain 
UV-prescription because E/SO galaxies contribute little to the total luminosity density at 
A < 3600A as long as redshifts of z < 1 are concerned in our analysis. 

In the galaxy evolution models, stars are assumed to be formed at a rate proportional 
to a power of mass fraction of gas in a galaxy such as C(t) = uf n (t). The rate coefficient 
v is chosen in a way to reproduce the present-day colors of each galaxy type assuming 
T G =15 Gyrs and the power index n is generally believed to lie in the range of 1-2 (see, 
e.g., Kennicutt 1989). We consider four different evolution models for spiral galaxies such 
as SI, S2, II, and 12 (for details see AYT). The symbol 'S' refers to the simple, closed-box 
model, while T to the infall model which allows for material infall into the disk region. 
The number attached to 'S' and T is the adopted value of power index n in C(t). The 
above four models are considered to cover a range of possible star formation histories and 
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are used to assess the uncertainty associated with the prescriptions of galaxy evolution in 
our analysis. 

In order to compare with the CFRS data, we set the normalization as 
£(T G ,\ B ) = 10 19 - 296 /i 5 o W Hz' 1 Mpc~ 3 , where h 50 = H /(50 km/s/Mpc), following 
Lilly et al. (1996). The contribution from each morphological type in C(Tq, X b ) is evaluated 
from the Schechter parameters (4> k , L* B k ,a k ) determined by Efstathiou, Ellis, & Peterson 
(1988) for the Anglo-Australian Redshift Survey (AARS), the Kirshner-Oemler-Schechter 
survey (KOS), the Revised Shapley-Ames sample (RSA), and the Center for Astrophysics 
redshift survey (CfA). The Schechter parameters, converted into the same 5-band system 
for the four redshift surveys, are summarized in Table 1 of Yoshii & Takahara (1989) 
and are used to calculate the proportion of local luminosity density among galaxy types 
as C k = 4> k L* k T(a k + 2). The results are tabulated in Table [TJ. When spiral galaxies are 
not classified into the four subclasses in the AARS, KOS, and RSA surveys, the relative 
proportion among the subclasses is assumed to be the same as that for the CfA survey. Fig. 
H] shows the luminosity density evolution as a function of time for each galaxy type assuming 
the CfA mix and Tg=15 Gyrs. In this figure the SI models are used for the evolution of 
spiral galaxies. Unless otherwise stated, the CfA mix is used below as a standard choice. 

3. Probing the Cosmological Parameters 

We consider three representative cosmological models: the Einstein-de Sitter (EdS) 
universe with (Qq, Ao) = (1,0), an open universe with (0.2,0), and a A-dominated, flat 
universe with (0.2, 0.8). The redshift of galaxy formation is assumed to be zp = 5. Different 
values of Zp hardly change our result as far as we restrict ourselves in z < 1 (see discussion 
in §|j). The Hubble constant h = f/o/(100 km/s/Mpc) is taken as 0.5, 0.6, and 0.7 for the 
three cosmological models in order to give a reasonable age of galaxies (12.1, 12.4, and 13.6 
Gyrs, respectively). Although the galaxy evolution models reproduce the present-day SED 
at Tg=15 Gyrs, we can use these models for the case of Tq= 10-15 Gyrs as well, because the 
luminosity evolution becomes very weak and changes the SED very little after t ~10 Gyrs. 

In order to compare with the CFRS data obtained for z < 1, we transform the 
theoretical predictions into those expressed in terms of z. In Fig. |2|, we show the comoving 
luminosity densities at 2800A, 4400A, and 1/xm as a function of z for the EdS universe 
(left panel), the open universe (middle panel), and the A-dominated universe (right panel). 
The CFRS data points are those from the "luminosity-function-estimated" values (Lilly et 
al. 1996) for 2800A (circles), 4400A (t riangles), and 1 /im (squares). The four theoretical 
curves for each waveband correspond to the variation of the evolution models of spiral 
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galaxies: SI (solid), S2 (dashed), II (dot-short-dashed), and 12 (dot-long-dashed). In this 
figure all the predicted and observed luminosity densities are shown after scaled to h = 0.5, 
and the theoretical curves are made to agree with the 4400 A data at z — 0. 

The behaviors of theoretical curves are not very different among the three cosmological 
models because the elapsed time between z = and 1 is roughly the same as 8 Gyrs. On 
the other hand, the observed luminosity densities, which have been derived from the direct 
observables of apparent luminosities and redshifts of galaxies, generally depend on the 
cosmological parameters. It is clear from Fig. |2| that the observed ^-evolution at z = — 1 
becomes flatter with decreasing Qq or increasing Ao, and eventually falls in agreement 
with the A-dominated universe. In sharp contrast, however, the observed ^-evolution is 
too steep to agree with the EdS or open universe, and such a discrepancy is much more 
considerable at 1/mi. We note that the /^-evolution is reliably modeled at 1/xm, offering 
an ideal cosmological test when compared to shorter wavelengths. This is because the 
luminosity density at 1/mi is always dominated by ellipticals and is quite insensitive to 
details of recent star formation history in spiral galaxies, as shown in the lower panel of 
Fig.l. 

For the purpose of quantitative comparison, we have calculated the slope index of the 
^-evolution, a = dlog£/dlog(l + z), averaged over the range of z = — 1 for each of the 
galaxy evolution models (SI— 12). We have repeated calculations with type mixes other 
than the CfA mix in Table [l], but the resulting a is changed at most by 0.08. Lilly et al. 
(1996) estimated the observed slope for the EdS cosmology as a = 3.90± 0.75(2800A), 
2.72 ± 0.5(4400A), and 2.11 ± 0.5 (1 fim). Using these observational errors which do not 
depend on Q and A , we have performed a x 2 analysis for a and found that the EdS 
and open cosmologies are inconsistent with the data with 99.86% CL. and 98.6% C.L., 
respectively. Assuming a flat universe (f2 + Ao = 1) with h = 0.6, the lower limit on Ao is 
obtained as 0.37 (99% CL.) or 0.53 (95% C.L.). Here, one of the galaxy evolution models 
(SI— 12) is used that gives the most conservative result in the x 2 analysis. 



4. Discussion 

Our conclusion supporting the A-dominated universe is drawn from using the AY and 
AYT models of galaxy evolution based on fairly standard assumptions on galaxy evolution 
and local properties of galaxies. However, some other scenarios of galaxy evolution have 
been proposed so far mainly to solve the so-called "excess" faint blue galaxies around 
B ~ 24. These include a scenario of widespread merging of galaxies (Broadhurst, Ellis, & 
Glazebrook 1992) and a scenario which invokes an exotic intermediate-redshift population 
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of dwarf galaxies which have ever faded away (Babul & Rees 1992). It is worth emphasizing 
that our analysis in this Letter is based on the luminosity density of galaxies but not on 
their number density. Luminosity density is generally less sensitive to mergers than number 
density, because increase of number of galaxies at high redshifts is compensated by the 
corresponding decrease of luminosity of individual galaxies. The exotic dwarf galaxies 
may significantly contribute in the number density, but it is doubtful that such a faint 
population occupies a large fraction in the CFRS galaxies at 0.5 < z < 1. In fact, their 
HST images exhibit the same range of morphological types as seen locally and the CFRS 
galaxies at z ~0.75 are broadly similar to the local galaxy population (Schade et al. 1995), 
indicating that exotic evolutionary processes do not dominate in luminous galaxies. 

The redshift of galaxy formation, though still uncertain, has been assumed to be 
zp = 5, because different choices give a convergent result below z ~ 1 and do not affect 
the cosmological test presented in this Letter. Recent progress in direct observations of 
star-forming galaxies at z > 2 however provides more useful constraint on zp. In Fig. 
|3], assuming the A-dominated universe with (f2 ,A ) = (0.2,0.8), we plot the observed 
star formation rates (SFRs) in the universe which are compiled by Madau et al. (1996), 
together with the predicted SFRs with zp = 5 and 7 assuming h = 0.7 (for details see 
Totani, Sato, & Yoshii 1996). The SFRs, observed or predicted, are shown after scaled to 
h = 0.5 in this figure. Thick lines show the SFRs from spiral galaxies, and thin lines from 
all types including elliptical galaxies. Since the absolute SFRs are uncertain by a factor 
of 2-3, we normalize the curves to agree with the data at z < 1. Inspection of this figure 
clearly indicates that zp is larger than at least ~ 4. According to Madau et al. (1996), the 
observed SFRs beyond z = 1 shown in Fig. should be considered as lower limits to the 
real values. However, since their analysis is based on a sample which covers a significant 
range of magnitude in the galaxy luminosity function, it may be safe to assume that these 
lower limits are not very different from the real values. If this is the case, the infall models 
of spiral galaxies with zp ~ 4-5 well explain the SFR trend at z ~ 2-4, whereas the model 
of elliptical galaxies overpredicts the SFRs unless they were formed at zp > 5. 

We point out that elliptical galaxies having high SFRs release metals out of the system 
into the intergalactic medium (IGM) through galactic winds. Miralda-Escude and Rees 
(1997) argued that the average metal abundance of Z ~ 2 x 10 -4 , which is observed in the 
Lya forests seen in quasar spectra, needs one Type II supernova per 5OOOM baryons in 
the IGM and these supernova progenitors should have emitted about 20 ionizing photons 
per each baryon, which are enough to reionize the universe. In our models, an elliptical 
galaxy with Mbaryon = 10 12 M Q produces about 2.5 xlO 10 Type II supernovae during the 
first 0.7 Gyr (Totani, Sato, & Yoshii 1996) and then releases metal-enriched gas of about 7 
% of Mbaryon into the IGM (AY). If the baryon density in the universe is taken as a typical 
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value of Vl b h 2 = 0.02 from the Big-Bang nucleosynthesis, the average metal abundance 
in the universe is given by Z ~ 1.7 x 10 4 /i 2 . Since this value agrees with the estimate 
by Miralda-Escude and Rees (1997), we can regard the elliptical galaxies as a source of 
supplying not only the metals in the IGM but also the ionizing photons which reionize the 
universe. 
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Time from Galaxy Formation [Gyrs] 



Fig. I. — Evolution of the comoving luminosity density at 2800A, 4400A and 1/im as a 
function of time for model galaxies of different morphological types. The galactic wind 
model from AY is used for E/SO (solid line), and the SI models from AYT for spiral 
galaxies of Sab (short-dashed), Sbc (long-dashed), Scd (dot-short-dashed), and Sdm (dot- 
long-dashed). The luminosity density at 15 Gyrs for the composite of all types is normalized 
to £(4400A) = 10 19 ' 296 /i 5 o W Hz -1 Mpc -3 , and the relative contribution among the types is 
adopted from the CfA mix in Table 1. 
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Fig. 2. — Evolution of the co moving luminosity density as a function of redshift for the 
composite of all galaxy types, shown for the Einstein-de Sitter universe with (fl , A ) = (1,0) 
(left panel), the open universe with (0.2, 0) (middle panel), and the A-dominated flat universe 
with (0.2, 0.8) (right panel). The CFRS data points (Lilly et al. 1996), scaled to h = 0.5, 
are shown by circles (2800A), triangles (4400A), and squares (1/zm). The four theoretical 
curves for each of 2800A, 4400A, and lfim correspond to the evolution models of SI (solid 
line), S2 (dashed), II (dot-short-dashed), and 12 (dot-long-dashed) for spiral galaxies. The 
galactic wind model for elliptical galaxies is used in common. The curves are normalized to 
coincide with the 4400 A data at z — 0. 
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Fig. 3. — Star formation rate (SFR) in the universe as a function of redshift. The A- 
dominated flat universe with (f2 , A ) = (0.2, 0.8) is assumed for both predicted and observed 
SFRs shown in this figure after scaled to h = 0.5. The data points are those compiled by 
Madau et al. (1996). The theoretical curves are based on the SI (solid) or II (dot-dashed) 
models for spiral galaxies and the galactic wind model for elliptical galaxies. Galaxies of all 
types are assumed to be formed at zf = 5 or 7, with h = 0.7. The SFRs from spiral galaxies 
only are shown by thick lines, while those from all types including elliptical galaxies by thin 
lines. The theoretical curves are normalized to agree with the data at z < 1. 



